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Abstract 

Presented are two kinds of integral solutions to the quantum Knizhnik-Zamolodchikov equations 
for the 2n-point correlation functions of the Heisenberg XYZ antiferromagnet. Our first integral 
solution can be obtained from those for the cyclic SOS model by using the vertex-face correspon- 
dence. By the construction, the sum with respect to the local height variables ko, k\, ■ ■ ■ , kin of 
the cyclic SOS model remains other than n-fold integral in the first solution. In order to perform 
those summations, we improve that to find the second integral solution of (r + l)n-fold integral for 
r £ Z>i, where r is a parameter of the XYZ model. Furthermore, we discuss the relations among 
our formula, Lashkevich-Pugai's formula and Shiraishi's one. 

PACS number: 02.30.1 



1 Introduction 

Thirty years have passed since Baxter's papers on the eight-vertex model Q and the XYZ spin chain || 
were published. In a series of papers || |], |) Baxter constructed the eigenvectors of the Hamiltonian of 
the XYZ spin chain by translating the eight-vertex model to the equivalent face model, the eight-vertex 
SOS model. 

Though many authors have tried to extract the results on correlation functions of the eight- vertex/XYZ 
model, very few results were obtained. The spontaneous polarization was conjectured by Baxter and 
Kelland ||, and the same results were reproduced by solving a set of difference equations, the quantum 
Knizhnik-Zamolodchikov equation (g-KZ equations) in Q . On the basis of the Z-invariance the integral 
formulae of correlation functions were obtained in the special case that the eight-vertex model can be 
factorized into two independent Ising model S . A free field representation of the type I and type II 
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JlO| vertex operators were constructed to express the correlation functions and the form factors of the 
eight- vertex/XYZ model in terms of those of the eight- vertex SOS model p~T| ] . 

There are several ways to address the problem of correlation functions of integrable models. The 
vertex operator approach p2| provides a powerful tool to calculate the correlation functions of the six- 
vertex/XXZ model the RSOS-type model p], p|, pfl. The bootstrap approach |l6| is based on 
the corner transfer matrix (CTM) method |ljj and enables us to derive a set of difference equations 
of correlation functions. Note that the CTM method can be applied to massive integrable model such 
as the XYZ Heisenberg chain in the antiferromagnetic regime. Kitanin et al jl8|, developed the 
algebraic Bethe Ansatz method to obtain domain wall correlation functions of the XXZ model, using 
the determinant structure of the partition function. 

In this paper we try to construct the correlation functions of the XYZ antiferromagnet by directly 
solving a set of difference equations, bootstrap equations. These equations are derived on the basis of 
CTM bootstrap in @. In a previous paper |]20f we did the same thing for both the bulk and boundary 
XXZ antiferromagnet. Let us cite some results concerning the bulk XXZ case from ]20| . 

The 2n-point correlation functions of the bulk XXZ antiferromagnet with the spectral parameters 
(C) = (Ci) " " i C2n) and the spin variables e = (si, ■ ■ ■ , £m) can be expressed up to a scalar function of 
C's in the following form: 

g ( ;\cy = n / d ™L ^({ Wa } aeA \oQ^({ Wa } aeA \cr, ^ = ±) (i.i) 



where 



A := {a\e a > 0, Ka^ 2n}. (1.2) 

The functions (u>|C) e is a meromorphic function defined below, and , I / i™' ) (u>|£) is a kernel with 
appropriate transformation properties |2(5, eq.(2. 32-33)] and recursion relations p(l eq.(2.34)]. The 



subscript a — ± in (1.1) specifies one of two vacuums of the XXZ model in the antiferromagnetic 



regime. The explicit expression of (w\() e is as follows: 



Q (n) Hcr = n & - w ^ n ^ - ^ n - ^r 1 , a.3) 

aeA \j=l j=a+l ) a,b€A 

\ ' a<b 

where Zj — £|. Here we notice the relation between our previous work [^(| and |22| , the latter 
in which the correlation functions of the massless XXZ model were obtained. Since the XXZ model 
has only one vacuum in the massless regime, the vacuum structure of the massless XXZ model is very 
different from that of the XXZ antiferromagnet. Nevertheless, the structures of the integral formulae 
for the correlation functions in the both regime are quite similar. In order to see such similarity, let us 
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(1.4) 



substitute x — e e , Q — e e/3j and w a = e 2eoia into (1.2). Then we have 

Q (n) HC) £ 

(o-l 2n 
JJshe(a -^-|) she(/%-a a - 
j'=l j=a+l 
2n 

x J]e" £ * (sheK-^ + l))- 1 . 

When |x| =1 (eg the expression ( |l.4| ) is equal to the meromorphic function Q n (a\f3) s in pl[ , 

the massless analogue of Q' n )(ii;|£) £ , up to a trivial factor. Furthermore, we should notice here that the 
corresponding expression for the eight-vertex SOS model which will be given in section 2 possesses also 
quite similar structure to (Oh. 

The XYZ Heisenberg antiferromagnet has two parameters q and p, where p — (—q) r . Lashkevich 
and Pugai || wrote down the n-fold integral formulae for the 2n-point correlation functions G in terms 
of those of SOS model F as follows: 

G(U1,--- ,U 2n \u Q ) = *fcJ("l-«0)® •••®*£r 1 ( U 2«- U °) ir ( U l'-" ,M2n|uo) fcofcr " fe2 "- (1.5) 

Here ui,--- ,W2n are additive spectral parameters, uo is an auxiliary parameter, and t^ ±1 (u) is the 
intertwining vectors introduced by Baxter || |J. The sum with respect to kj (1 ^ j ^ 2n) should be 
taken over fcj = fcj-i ± 1. Concerning the sum with respect to ko, that depends on the value of r. In 
(L5) r > 1 is supposed. When r is a rational number of the form 2r = N/N' (N, N' coprime), we 
should take the sum over Z/7VZ, otherwise over Z. 

On the other hand, Shiraishi |2^] constructed the formulae of the correlation functions of the 
XYZ model without using the vertex-face correspondence. There are three main differences between 
Lashkevich-Pugai's formulae and Shiraishi's ones. Shiraishi restricted himself to the case r = |, while 
r > 1 in (1.5). Furthermore, the 2n-point correlation functions a la Shiraishi are of the 2n-fold integral 
form, while n-fold in the former case. These two differences seem disadvantage of Shiraishi's formulae 
in a sense, however, the third difference gives advantage. Shiraishi constructed the type II and type 
I vertex operators of the eight-vertex model as intertwiners of the g-deformed Virasoro algebra |f24j| , 
without transforming those of SOS model via vertex-face correspondence. Thus, the sum with respect 
to kj in ( jl.5| ) are not needed. 

In the present paper we wish to synthesize the advantages of the two formulae. Namely, we construct 
(r + l)n-fold integral formulae for the 2n-point correlation of the XYZ model with r G Z>i, by using 
the vertex-face correspondence, and performing the sum with respect to kj. 

The rest of the present paper is organized as follows. In section 2 we formulate the XYZ antiferro- 
magnet and the corresponding cyclic SOS model. In section 3 we present the bootstrap equations and 
construct integral formulae for correlation functions for the cyclic SOS model. In section 4 we further 
present integral formulae for correlation functions of the XYZ antiferromagnet. In section 5 we give 
some concluding remarks. In Appendix A we list some properties of the i?-matrix of the eight-vertex 
model. 
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2 The XYZ antiferromagnet and the cyclic SOS model 



2.1 The Hamiltonian and the i?-matrix 

The Jacobi theta functions with the characteristics a, b € R are defined by 



a 
b 



(u; t) := ^ exp {ir\^l(m + a) [(to + a)r + 2(u + b)] } 
The ith theta functions are defined as follows: 

(u;t), 2 (w;t)=^ 

(u;t), 6> 4 (u;t) = i? 



(2.1) 



9i(w;r) = i? 



1/2 
-1/2 



1/2 





1/2 



(u;t). 



These theta functions have the infinite product forms, e.g. 



l(w;t) 



-lq 4 e 



^(e 2 ^ 1 ™), 9 2 (u-t) = qZe-^^Q^i-e 2 ^ 1 ™), 



where q = exp(V~ Ittt) and we used the standard notation 

Q p (z) := (z;p) 00 (pz~ 1 ;p) 00 (p,p) 00 , (a;pi, • • • ,p n ) c 



fcl 
aPi 



■P k n)- 



ki>0 



In this section we consider the XYZ spin chain in an infinite lattice H [jjj 



(2.2) 



Here <tJ, txj and denote the standard Pauli matrices acting on j-th site, and we restrict ourselves to 



the antiferromagnetic regime |T| < 1 and A < — 1, where 



r = 



1-7 



1+7' 
A = - 



7 = 



Q2( V — le , 2y^ler \ 



> 0, 



1 9i(o ; ^^)(e 2 e 3 )(^; 2 -^) 



1 + 7 (^ 3 )(0; 2^i2)02(^ ; 2^12) 

The condition |T| < 1 and A < — 1 is equivalent to the one such that e > and r > 1. The antiferro- 
magnetic XXZ Hamiltonian can be obtained by taking the limit r — > 00. In this limit we have 7^0 
and A — > — (x + x _1 )/2, where < x = e _e < 1. As well as known, thus obtained XXZ Hamiltonian 
commutes with the quantum affine algebra U- X (slz) [ p~2| . Let = Cv + + Cu_ be a vector represen- 
tation of f7_ x (s[ 2 ). Then the limiting XXZ Hamiltonian formally acts on V®°° — ■■■^V^V®---. 
In [^2| the space of states I/® 00 was identified with the tensor product of level 1 highest and level — 1 
lowest representations of U-xish). It is very likely that the same structure survives in the nonlimiting 
generic case. In other words, the XYZ antiferromagnet is expected to have the symmetry described by 
the elliptic affine algebra A q , p (s\.2) JUi Uli where q = —x and p = x 2r . 
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The XYZ Hamiltonian ( |2.2| ) can be obtained from the transfer matrix for the eight-vertex model, 
by taking logarithmic derivative with respect to the spectral parameter £. The R- matrix Rgv(C) € 
End(C 2 £g> C 2 ) of the eight- vertex model is given as follows: 

#8v(C) = w±{u)I ® I + w 2 {u)a x ® a x + w z {u)v v ® a v + wi(u)<j z ® <r z , 
where z = ( 2 = x 2u and 

a v 2r r ' er ' 



w a {u) 



2r ' er 



In what follows we denote Rsv(C) an d -fiW(C) by -R(C) an d R(()> respectively, when there is no fear of 



confusion. In the above expression we modify Sklyanin's paramctrization |2?H such that (2^2) coincides 
with Baxter's one |r^| in the principal regime. The normalization factor 

(2.4) 



«(0 
s(C) 



[i] 



1 ( ' 1* ' T* T 1 ^ ( IP " iy "T* ^ ^ T* y* T* T* ^ ^ T 1 ~ ^* "I - 2 J- • /t>4 ™2r \ 



f™4 ? -l. ™4 ~2rA /^2 r .^4 ™2r^ (™2r r 



2r~-l. ~4 ~2r 



5 (a; 



2r+2 r . ™4 ^.2r s | : 

,6, X , ^ joo 



is chosen such that the partition function per site is unity. In other words, the factor is the partition 
function per site of the unnormalized model defined by R(C)- For later convenience we also introduce 
the symbol {it} by 

{ u } = x^- u e x ,r(-x 2u ). 

Note that 

^) = /fcxp (-f ) [u], 9^- ^§1) = exp (-f ) {u}. 
Let C 2 = Cv+ © Cv_ and introduce the matrix elements of the i?-matrix as follows: 



R(()v £l ®v e2 = ]T v £li ®v £ ,R(()ii = v £ ,®v e . 2 R({)£i 

K{ ^> e'e>=± 



(2.5) 



e' lt s' 2 =± 

Then the nonzero entries are given by 

R{QXX = R(0~- = o(C) 



R(C)X- = R(0-+ = KO 



(W(o ; ^)(M 4 )( 



2r 1 er 



R(0±+ = R(()+± = c(C) 



(W(0;^)(<M4)(£;^) ^ 
^)(W(^;^) 



(2.6) 



(W( 



2r 5 er 



i?(Otl = i?(o;; = d(c) = - 



(W(O;^)(0 2 3 )(£;^) 



2r 



Assume that the parameters it, e and r lie in the so-called principal regime |^7| : 

e > 0, r > 1, - 1 < u < 0; < p < x < C -1 < 1, 



where x = e e , p = x and £ = a;". The main properties of the i?-matrix are listed in Appendix A. 



G 



Yas-Hiro Quano 



2.2 Vertex- face correspondence 



In order to construct the eigenvectors of the eight-vertex model, Baxter introduced the following inter- 
twining vector j|, |4|] : 



t k k±1 (u) = f(u)t k k±1 (u), 

t k k±1 (u) = (^T) fe (M 4 )fe^ 



v+ + (v/=T) 



-kt 



no 3) 



2r 



The normalization factor f(u) satisfies the relation 

C 2 [u]f(u)f(u-l) = l, 
The explicit expression of f(u) is as follows: 



C=^e 

7T 



(x 2r ,x 2r ) 
(-x 2r ,x 2r ) 



X 2r + ' 2 r U + 4, (x* 4 Z ] X^ , X^ ) qo {x 2V + 2 Z J , ) ^ 



(2.7) 



(2.8) 



(2.9) 



CV(a; 2r ;^ 2r )oo (x 2 z;x 4 ,x 2r ) 00 (x 2r z;x 4 ,x 2r ) c 

The SOS model is a face model [|| which is defined on the square lattice with a site variable ki e Z 
attached to each site i. We call ki local state or height and impose the condition that heights of adjoining 
sites differ by one. Local Boltzmann weight are assumed to be the function of the spectral parameter 

c d 



u and to be denoted by W 



u , which is given for a state configuration 



round a 



face. Here the four states a, b, c and d are ordered clockwise from the SE corner. The Boltzmann weights 
of this model are assumed to be the function of the spectral parameter u and the nonzero Boltzmann 
weights are given as follows: 



W 



w 



k±2 k±l 
k±l k 



1 



k(u) ' 



k±l 



k±l 



k±l 



kT 1 



1 [#T«} 
k(u) [u + l]{k}' 



(2.10) 



u \ 



1 [u]{fc±l} 

k(u) [u + l]{k}' 



The normalization factor k(u) given by ( |2.4| ) is chosen such that the partition function per cite is unity. 
Then we have the so-called vertex- face correspondence H: 



R(ui - u 2 )t c b {ux) (8 t b a {u 2 ) =Y. W 



ui-u 2 \ t d a {ui)®t c d {u 2 ). 



(2.11) 



The RHS's of the last two equations in ( 2.10| ) contain the factors like {k}, while the corresponding 
Boltzmann weights of the ABF model |2£| are expressed in terms of the factors like [k] . Suppose that 
r is a positive integer greater than 1. Then the ABF model is called a restricted SOS (RSOS) model, 
since [k] = for k £ rZ and therefore the local states can be restricted as fc = l, — 1. On the other 
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hand, the present face model (2.10) is a cyclic SOS model because {k + r} = {k}. In || the RSOS type 
weights were used so that the following regularization was required. In order to avoid the pole resulting 
from [k] in the denominator, k £ Z + S should be assumed with some real S, and the limit 5 — ► should 
be taken after all calculation M. In our case we need not such regularization because {k} ^ for k e Z. 



We also notice that the intertwining vector (2.7) is closely connected with 2r dimensional cyclic 
representation of Sklyanin algebra Actually, the L-operator defined via R^yLL = LLR S v 



relation can be factorized into the intertwining vector (2.7) and its dual vector 



3 Bootstrap equations in the cyclic SOS model 
3.1 Integral formulae for the cyclic SOS case 

As preliminary, let us consider the bootstrap equations and correlation functions for the cyclic SOS 
model. Let 



(3.1) 



for \kj — kj—i\ = 1 (1 < j < 2n) and ko = ki n = k, be the 2n-point correlation function of the 
eight-vertex SOS model. Then these functions satisfy the following three CTM bootstrap equations 
E3: 



(" ' ' > 0+1)0) ' ' ' )' 

fe' \ k j k i-i 



' kj i_ kj kj _|_ 2_ * 



o/o+i ,o,o+i.---r fci - ifc ^ +i - 



(3.2) 



F^(Ci,--- X2n-i,x 2 ( 2n ) k --- k ^ k = ajj^F^iC^Ci, - ■ ■ ,C2„-i) fe2 "- lfe - fe2 "- 1 . (3.3) 

Fi n) (Cl,"" ,(2n-2,(2n-l,X- 1 C 2n . 1 ) k - k ^ k ^ k = 5 ^± F ^ (0 , ■ • • , ^n-^ (3.4) 



r<y " ! S.2«-2j S,2«-l) — x (,2n-l) — O, 



k 2n -2,k v /_ 1 | fc | -tV lU, - • • > <,2n-2) 



Set 



( f~\ ^1 ' ' ' &2n — 1 k 



F (n) {C) kk 1 -k 2n - 1 k = Cn Cp^g( Zj /z k )xF™{Q 
l^j<k^2n 

Here c n is a constant which will be determined below, and the function g(z) has the properties 



g(z) = g(x-*z- L ), k(C) = C 



(3.5) 



(3.6) 



(3.7) 



The explicit form of g(z) is as follows: 



^6„. T 4 4 2r\ / 2 -1 . 4 4 2r\ („,2r+6 _. ™4 4 2r\ / 2r+2 -1. 4 4 2r\ 

j oo \+ ~ , a- , a/ , ct/ j oo V ) ) ) *^ / oo \ , x , x , x 700 /q q \ 

^v 2 ' — f T 8,. -r-4 ™4 T ,2r N i TZIZ-L. Zi ™4 ™2r"\ /™2r+4 r . ™4 ^4 ™2r'^ ^™2r r -l. ~4 T 4 „,2r^ ' 
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( n ) 

In order to present our integral formulae for F a (£) let us prepare some notation. Let 

A := {a\k a = k a _ 1 + 1, 1 sC a < 2n}. (3.9) 

Then the number of elements of A is equal to n, because we now set &2n — k — k. We often use the 
abbreviations (w) = (w ai ,--- ,w an ), (w') = {w ai ,--- ,w an _ x ) and (w") = {w ai ,--- ,w an _ 2 ) for aj e A 
such that a\ < ■ ■ ■ < a n . Let us define the following meromorphic function 



Q(")( W \0 kkl - k2n - lk = l\{Va-U a + h 2 -k a } [HlVa-Uj-ft J] 1 



a— 1 2n 



2n-l 



(3.10) 



a<6 



where w a = x 2Va . Here we should notice that the structure of the expression ( 3.1C| ) is quite similar to 

We wish to find integral formulae of the form 

j?W _ -Q ^ ^ ^{w\QQ^{w\C) kkl "' k2n - xk . (3.11) 

Here, the kernel has the form 



aeA JC a 2nV-lw a 



2/i 



^Ho=<" ) Honn^^^^(^) n (3.12) 

a£Aj = l 3 l^j<kf^2n 



where 



/„2r+3 „,4 „2r\ I „2r+3 7 -l . „4 „2r\ 



f T -.»4 1 ,2r 1 i C_--l._4 -Sri ' (3.13) 



For the function we assume that 

• it is anti-symmetric and holomorphic in the w a € C\{0}, 

• it is symmetric and meromorphic in the Q G C\{0}, 

• it has the two transformation properties 



^MC^&O = ax -^ + ^Y[ Wa flC-\ (3.14) 



■d { J l) (w',x i Wa n \0 _ _4„TT z . 



a£A j=l 



it satisfies the following recursion relation 

2n-2 _ t (3.16) 



a; 



n— (3-i5) 
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Here we also fix the constant c n as follows: 



{x 2 r O0 {x e r 0O {x^+ 2 r O0 w r+6 } n 



(3.17) 



where 



{^■}oo • (^"5 X 7 X j X ) 



The function (HO i s otherwise arbitrary, and the choice of d^ 1 ' (w\Q corresponds to that of solutions 
to (3.2-3.5). The transformation property of (w|C) implies 



a(™)/ 



^\w',x 4 w a JC) 



2n 



n 



n 



U 2 r 



U2r, 



(3.18) 

(3.19) 



The integrand may have poles at 



±(i+ik+2ri) z . (i ^ j ^ 2n, k, I e Z >0 ), 
(6 < a), 
(6 > a). 



X 2 Wb 



(3.20) 



a; 2 Wf, 



We choose the integration contour C a with respect to w a (a £ A) such that C a is along a simple closed 
curve oriented anti-clockwise, and encircles the points x 1+Ak+2rl Zj(l ^ j ^ 2n,k,l £ 2^o) an d x 2 Wb 



(b < a) but not x 



-l-ik-2rl 



z j(l ^ j ^ 2n, fc, Z € 2^o) n or a; 2 ui(, (6 > a). Thus the contour C a actually 



depends on z^s besides a, so that it should be denoted by C a (z) = C a {z\, • ■ • , Z2n), precisely. The LHS 



of (3^3) refers to the analytic continuation with respect to C,2n- Nevertheless, once we restrict ourselves 



to the principal regime < p < x < < 1, we can tune all C a 's to be the common integration contour 
C : \w a \ = x^ 1 because of the inequality xzj < x" 1 < x~ 1 Zj. 

We are now in a position to state the following proposition regarding the correlation functions and 
the bootstrap equations in the cyclic SOS model case: 



Proposition 1 Assume the properties of the function below ( 3. IS ) and the integration contour 



C a below i p.2(\ ). Then the integral formulae \3.6j , \3.l\ ) with Sj , \3.1<\ , \3.13j , \3.13j ) solves the three 



equations (3.i-3.l) 



The proof of Proposition [T] will be given in the subsequent subsections. 



3.2 Proof of the M^-symmetry 



Let us first prove (3.2). For that purpose we have to consider four cases according as kj — fcj-i = ±1 
and fcj+i — kj = ±1. 



Suppose that kj — fcj-i = fcj+i — kj = —1. Then the relation (3.2) holds, because the integrand 



(„) 

F a (£) is evidently symmetric with respect to Q and Cj+i- 
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Let (fcj-i, kj, fcj+i) = (to, m — 1, m) for some to. Then the relation ( |3 .2| ) reduces to 

h 0+1 : 

i+i) ' ' ' 



■Fff • " ) 0+1 ! Cj! 



[l]{m + Mj -M j+ i} -(n) ...,„,„ . ,,„. 

[n i - %+1 +l]{m} CT ^••^'■^ 
[uj-u ]+1 \{m + l} —(n) f , , .„„„+,.„,. 



(3.21) 



Note that the set of the integration variables in the second term of the RHS is different from the other 
terms. Since w a = x 2Va 's are integration variables, we can replace both Vj and in the integrand by 



v. After that, the relation ( 3.21 ) follows from the equality of the integrands. In this step we use 



r i ii [1]{to+ u,- — r -i , r 

{v — Uj + I — m\[v — Uj+i — ij\ = 7 -tj — t-\ v " u j+i + 2 — ~ 



itj — + l]{m} 
Uj — Uj+i]{m — 1} 



{v — Uj — ~ — TO,}[Hj + l — 



J 2 J 



[uj-u J+ i + l]{TO} 

Suppose that (fcj-i, fcj, fcj+i) = (to, to + 1, to) for some to. This case can be proved in a similar way 
to the previous case. Here we use 

[l]{m-u j +u j+1 } 



{v - uj+i - | - m}[uj -v - §] 



% — Uj+i + 1]{to} 
Uj — Uj + i]{to + 1} 



[ Uj -u j+1 + 1]{to} 



{v — Uj — \ — m}[uj + i — v — |] 
{w — Uj +1 + i — m}[v — Uj - |] 



Finally, let (%_i, fcj, %+i) = (to — 1, to, to + 1) for some to. Then the integrand of the RHS of (3.2) 
contains the factor 

{wj-Uj + 5-m}[tt i+ i-i;j-|]{u i+ i-u i+ i-|-m}[uj + i-Uj-|] 
7(uj,Uj + i;Uj,Uj + i) = - — : . 

[Vj - Vj+l + 1J 

The corresponding factor in the LHS should be equal to /(itj+i, Uj] Vj,Vj+{). Thus the difference between 
both sides contains the factor 

I(uj,u j+ i;Vj,Vj + i) - I(uj + i,Uj-,Vj,Vj + i) = {m}[uj - u j+ i]{vj + v j+ i - u 3 - u J+ i - to}, 

which is symmetric with respect to Wj = x 2Vj and u>j+i = x 2v ^ +1 . Since ■ (iy|£) is antisymmetric 



with respect to u> a 's, the relation (3.2) in this case does holds 



3.3 Proof of the cyclicity 

In the proof of the cyclicity (|3.3|) we have to consider the two cases fen-i = k ± 1. First let fen-i = 



k + 1. When the integral ( 3.11 ) is analytically continued from (2n — x U2n to x 2 C,2n — x U2n+2 , the 



points x 



l+4k+2rl 



Zm and x 



-l-4k-2rl 



%2n (k-l S Z^o) move to the points x 5+ + Z2 n and 



,3-4k-2rl 



Z2m 



respectively. In the LHS of (3.3) the point w a = x Z2 n {v a — U2« + 4) are therefore outside the integral 



contour C' a = C a (z' , x 4 Z2n)- Nevertheless, we can deform C' a to the original one C a — C a (z) without 
crossing any poles. That is because the factor 



W [U2n ~ Va + §], 



aeA 



contained in Q^ n \w\(^' , x 2 ^2n) k ' k+lk cancels the singularity at w a = x 3 Z2 n - Thus the integral contours 



for both sides of (3.3) coincide 
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Furthermore, by using (3.1 8|) we obtain 



^Hw\c,x 2 c 2n )Y[[v a - U2n -i]^a¥j i Hw\oY[[u 2n -v a + ^ n (r-)' > 

a£A aeA j=l ^ ^ 2n ' 



which implies that the integrands of both sides of (3.3) coincide and therefore the relation (3.3) holds 
when k 2n -\ = k + 1. 

Next let fc2n-i = k — 1. In this case we make the rescale of variable w 2n x 4 w 2n (v 2n i— » v 2n + 2) in 



the LHS of (3.3). Then the integral contour with respect to w a (a £ A\{2n}) will be C' a = C a (z', x 4 z 2n ), 
and the other one will be C = C 2n (x~ 4 z' , z 2n ). For a S A\{2n}, we can deform the contour C' a to the 
original C a without crossing any poles, as the same reason as the previous case. The integral contour C 
encircles x~' i+4k+2rl Zj and x 1+4k+2rL z 2n , but not x~ 5 ~ 4k ~ 2rl zj nor x~ 1 ~ 4k ~ 2rl z 2n , where 1 < J < 2n— 1, 
k,l e Z^q. Since Q (n) (w\C, x 2 C 2n ) k '" k ~ lk contains the factor 



2n ~ 1 [u — v + -] 

J(w'x 4 W 2n \C ,X 2 ( 2n ) = {«2n - U 2n + \ - k} [V2n ~ Uj + §] ]~J " 2 



J'=l 



fo - v 2r. 



(3.22) 



the pole at w 2n — x 3 zj (1 ^ j ' ^ 2n— 1) disappears. Thus we can deform the contour C to the original 



one C 2n = C' 2n (z) without crossing any poles. Thus the integral contours in both sides of (3.3) coincide. 



Replace the integral variables such that (w' ,w 2n ) i— > (w 2n ,w r ) in the RHS of (3.3), and compare the 
integrands of both sides. Note that ( 3.22 ) describes allw 2n - and Z2n-dependence of Q( n \w\£', x 2 Q 2n ) k '" k ~ lk . 
From (3.18, 3.19) and the antisymmetric property of \E'^ with respect to w a 's, we have 



2n-l 



^J l Hw',x 4 w 2n \C,x 2 ( 2n )J(w , x 4 w 2n \C,x 2 (; 2n ) = a^\w 2n ,w'\C)J(w 2n ,w'\C) JJ , 



V 



, (3.23) 



which implies that the integrands of both sides of (3.3) coincide and therefore the relation (3.3) holds 
when fen— l = k — 1. 



3.4 Proof of the normalization condition 



Let us prove (3.4) and (3.5). The factor g{z 2n -\/ z 2n ) has a zero at (, 2n = ±x £an— l- On the other 
hand, the two points xz 2n and x z 2n —i are required to locate opposite sides of the integral contour 
C a , so that a pinching may occur as C, 2n — > ±x _1 C2n-i- If there is no pinching the correlation function 
j#0 (£)*•••* vanishes at ( 2n = ±x~\ 2n - 1 . 

Suppose that (k 2n - 2 , k 2n -i, k 2n ) = (fc+2, k+1, k). In this case the factor [v a — 1*271—1 + 5] contained in 
Q( n ) (i D \Q k --- k + 2k + lk cancel the poles of i}>{w a / z 2n -\) at w a = x~ 1 Z2 n -i, and therefore no pinching occurs 



as C, 2n — > ±a; -1 £2n-i- Thus the conditions (3.4) and ( |3.5| ) hold when {k 2n - 2l k 2n -\, k 2n ) = (fc+2, fc+1, fc). 

When (k 2n - 2 , k 2n _i, k 2n ) — (fc — 2, fc — l,fc), no pinching occurs for the integral contour for C a 
(a 6 A\{2n — l,2n}), as the same reason as the previous case. Concerning the integral with respect 
to w 2n -i and w 2n , there is no singularities at w 2n -x — xz 2n and w 2n — xz 2n and consequently no 
pinching occurs actually. In order to see such vanishing singularities, we first notice that the function 
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Q ( - n \w\() k "' k 2k lk contains the factor: 
q(v 2 

{t>2n-l ~ M2«-l + | ~ fe}[M2n ~ V2n-1 ~ \]{v2n ~ U 2n + \ ~ k}[v 2n ~ U 2n -1 ~ \] 

[V 2n -i - V 2n + 1] 

Because of the antisymmetric property of with respect to u> a 's, the zero of 'd^'iw", w 2n ^i, xz 2n \() 
at w 2n -i — xz 2n cancels the poles of ip(w 2n -i/ z 2n ) at the same point. Furthermore, the poles of 
^(w 2n j ' z 2n ) and ip(w 2n / z 2n _i) at w 2n — xz 2n and z 2n = x~ 2 z 2n ^\ are canceled by the zeros of 
g{z 2n -\/ z 2n ) and (w" , w 2n -i, xz 2n \() k ••• fe - 2fc - lfc a t the same points. The latter cancellation can 
be shown as follows. Note that the integral is invariant as we replace (w" ,w 2n _i, xz 2n \() k "' k ~ 2k ~ lk 
by its antisymmetric part with respect to w 2n ~i and xz 2n . Accordingly, we can replace the factor 
q(v 2n -i,u 2n + ^;U2n-i,U2n) by the following factor: 



q(v 2n -i,u 2n + \;u 2n -i,u 2n ) - q(u 2n + ^,v 2n ^i;u 2n ^i,u 2n ) 

{u2„-u 2 „-i+2-fc}[l 

U9„ I H 

[«27l+n-f2l, 



, ;; , , , ., , , r , {«2 n -«2T 1 -l+2-fc}[l]{l>2 n -l-U2r 1 + ^-*:}[u2, 1 -l-«2 re -l-^] 

{v 2n ~i - u 2n -i + 2 - k}{l - k\[u 2n -\ - u 2n \ H -g = 



that vanishes when z 2n = x 2 z 2n -\ (u 2n — u 2n ~\ — 1 or u 2n — u 2n -\ — 1 — Thus there is no 

singularity at w 2n — xz 2n as C2n — * ±2; C 2 n—i- The same thing at w 2n -\ — xz 2n can be easily shown in 



the same way. Hence the conditions (3^) and (3J5) are verified when (k 2n -i, k 2n -i, k 2n ) = (k— 2, k— 1, k). 

Next let (k 2n - 2 , k 2n -i, k 2n ) = {k,k— 1, k) and consider the limit (^ 2n — > x _1 C2n-i- In this case there 
is no pinching for the integrals with respect to w a (a S A\{2n}). Let C denote the integral contour 
with respect to w 2n such that C encircles the same points as C 2 n does but xz 2n . Note that no pinching 
occurs with respect to the integral along C because both the two points x~ l z 2n _\ and xz 2n lie outside 
the contour C. Thus the integral with respect to w 2n along the contour C 2n can be replaced by the 
residue at w 2n = xz 2n . 

In order to evaluate the residue, the following formulae are useful: 



r , z 2n -i \ z' xz 2r 



z 2n ^x Z2n-l \ Z 2n J \Z2n-\ , 

1 \X ^\qo\X }oo{*£ }oo{^ ~^ }oO 

(-.2. „2r\ („.2r-2. „,2r\ r T> 21 7^61 f T ,2r+21 J^2r+6\~ 



(3.24) 



P dw 2n ; ( w 2n \ 1 



g(z)g(x 2 z)^(xz)[u + 1] = x'^— ^ u+1 \x 2r ; x 2r ) oc . (3.26) 



*(!)* 



x 



-(«-t,-l/2)( a .2r ;iB 2r) c 



[u-U-i] = 5w / 2^ _ " ( 3 ' 27 ) 

(a; w/ z ; x l ) ^ [xzj w ; x z ) ^ 
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Using these, the condition (3.4) with (fc2n-2, fen-ij k 2n ) — (k,k — 1, fc) reduces to 



2r+4\ 2n ~ 2 



j'ooj^ j'oo'l*^ }oo{-£ } 
/ T 21 J" -61 f r 2r+2| f T 2r+61 



nc~n 



c D a -«2n-l + 4 

* Q x 2(xw a /z 2 n — l I 



which is valid under the assumption of (3.16) and (3.17). The relation ( |3.5| ) with (k 2n ^ 2 , k2 n —i, k 2n ) = 
(A;, fc + 1, fc) can be similarly proved. 

When (k 2n -2, fc 2n -ij ^2n) = (k, fc+1, fc), the only difference from the previous case is that the rational 
function (w\Q'" kk+lk contains the factor 



{V2n-1 - U2 n -1 — \ — k}[u 2n - V 2n -1 - \] 



{u 2n - U 2n -1 ~ k}[l] 



•W2n-l=XZ2n 



{fc+l} 

in the present case, while the corresponding factor in the previous case is 
{v 2n -u 2n + \~ k}[v 2n - u 2n _ 1 - i] 



{fc + 1} 



(3.28) 



{k-1} 



[u 2n - u 2 „_i] 



(3.29) 



Since ( 3.28 ) is equal to ( 3.29 ) as C, 2n — x 1 Q 2n -i, the condition (3A) with (k 2n - 2 , fen-i, k 2n ) = (fc, k 



1, fc) follows from the previous case. The condition (3.5) with (k 2n ~ 2l k 2n -\, k 2n ) = (fc, fc + 1, fc) can be 



similarly proved. After all, the conditions (|3.4|) and (3.5) were componentwisely proved 



3.5 Non-trivial theta function 



In subsections 3.2-3.4 we proved Proposition [j]. The key point in proving (3.2) was the 14^-matrix 



symmetry of the rational function Q( n )(u>|£) fe '" fe . We observed that the other two conditions (3.2) 



and (3.4) hold under the assumption of the transformation properties and the recursion relation of 







(n) 



Actually, you can easily find that (3.14), (3.15) and (3.16) are sufficient conditions of (3.3) for 



k 2n -\ = fc + 1, (3J5) for fen-i = k — 1 and (3.4), respectively. 

In this way we obtained the integral solutions to CTM bootstrap equations for correlation functions 
of the cyclic SOS model, with the freedom of the choice of da . Now we wish to present an example of 



■ff^ satisfying all the properties given below Q3.13p 



2 ii 



tf^Mc) = e aa -ax n w- 1 no n $ 



-(n-D(l-i) 



a£A 



n 

a<b 



x Vb - v «e X 2(w a / Wb ). 



(3.30) 



You can easily see that (|3.30| ) satisfies all the properties of symmetry with respect to Q's, antisymmetry 



with respect to u> a 's, ( |3.14 ), ( 3.1 5| ) and ( |3.16 ). We also notice that there is actually no poles at w a = x' 



Wf, (b > a) when we fix d„ to (3.30) because of the factor Q x 2(w a /wb 



(b < a) and w a = 

Furthermore, we should like to notice that this example of the function \w\() is quite similar to the 
anti- ferromagnetic XXZ analogue [^0| . 
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4 Integral formulae for the XYZ antiferomagnet 

4.1 Correlation functions and difference equations 

Let us introduce the V® 2n -valued correlation functions 

4 n) (Ci»'--.C2»)= ^v et ®--.®v e ^GW(( lr -.,& n ) ei '" l » n , (4-1) 
where a = ± signifies one of two vacuums of the XYZ Heisenberg antiferromagnet. Let 

= Efl, i ---E^ (4-2) 



where E~?) , is the matrix unit on the j'-th site. Then the correlation function (J4 J 1|) gives the expectation 



value of the local operator (4.2) by specializing the spectral parameters as follows: 



(0) CT = GW(ar 'C,--- ,x~\Xr-- ,C)- £ »-- £ ^- £ ». (4.3) 

In what follows we often use the abbreviations: (() = (0, • • • > (2n)> (CO = (0, ' ' ' , Can-i), (C") = 
(Ci 3 • * • iC2n-2), (s) = («i,--*-?2n) 5 ( z ') = (zi, • • • Z2n-i); and (e) = (ei.--e 2n ), (e') = (ei • • • e 2n -i), 
(e") = (ei • • • £271-2)- O n the basis of the CTM (corner transfer matrix) bootstrap approach, the 
correlation functions satisfy the following three conditions Q: 

1. R-matrix symmetry 

P„ +1 G^(... ,0+1,0,- ••) =R jj+ i(C j /Cj+i)G^(--- ,0,0+1, •••) (l<j<2n-l), (4.4) 

where P(x ® y) = y ® x. 

2. Cyclicity 

Pl2 ■ ■■P2n-12nG^((>,X 2 ( 2n ) = aG^^O- ( 4 -5) 

5. Normalization 

Gi" ) (c",c 2 r l -i,c 2 „)ic 2 „^- 1 c 2 „_ 1 - Gir^co®^ (* = ±), (4-6) 

where u s = v + ® V- + sw_ <8> f+- 

These three conditions can be componentwisely recast as follows: 



ci n) (--- ,0+1,0,- ••r™- 
2 i?(o/o + i):o:;^ n) o--,o,o+i,--: 



(4.7) 



G^(C',x 2 C 2 „) £ ' e2 " = aG^(C 2 „,C') £2 " £ '. (4.8) 
Gi" ) (C",C2n-i,^- 1 C 2n -i) £ " £ — 1£2 " = egrf /a ^»- 1 -Hr a »,o^" 1) (C /, ) e " (* = ±)- ( 4J ) 
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Combining (4.7) and ( [4.9[ ) we obtain another expression of the normalization condition: 

^ £ (i- S )/2 G w (c // iC2n _ iiS;rC2n _ i)e " e - e = Gi^gy, 

where we also use 



Risx' 1 ) = 



( \ 
OslO 
1 s 
y J 

Note that the first two conditions imply the difference equation of the quantum KZ type pit] of level 



(4.10) 



(4.11) 



where Tj is the shift operator such that 

T i F(C) = J , (Ci,--- ,x~%,>-- .Can), 
for any 2n-variate function i* 1 . When the definition of Tj is replaced by 

T j F(Q=F(h,--- ,i'+ 2 (i,- ,Ca»). 

and all the arguments x~ 2 Q/(k in the first line of the RHS of ( 1.10] ) are also replaced by x l+2 Q/(k, the 
difference equation ( 4.10 ) is called the quantum KZ equations of level I. 

Remark. As a result of Z 2 -symmetry of i?-matrix, there are two ground states in the XYZ anti- 
ferromagnet. Let us specify the two ground states by i = 0, 1, and denote the correlation function on 
the i-th ground state by G^ (£) . Representation theoretical speaking, the correlation function G\ (£) 
refers to the trace of type I vertex operators on the irreducible highest weig ht module of A g ,p(sl 2 ) 
|§. The CTM bootstrap approach suggest that both G^\() and G^ n) (C) will appear in the cyclicity 
condition as follows: 



P 



2n-12 



nG\ U \Cl, ' ' ' , Qn-l, X 2 (^2n) — G\J i (C2n, Cl) ' ' " ! C2n-l) 



{n) 



Thus we introduce G^ (£) — Gg (£) + crG^ such that the second equation (4.5) involves only 
G^(C). 



Thanks to (3.7) the first two equations (4.7-4.8) are rephrased in terms of G^ (£) and R(() as 
follows 



g { ;\--- ,0+1,0, •••r £i+i£ '-- = J2 Riti/MTJ+y? (••• ,G,o+i,--o e ^ +i , (4-12) 

^— * j ' j+i 



G| n) (C,* 2 c 2 „) £ 



= aGi n) (C 2 „,C') £ 



0=1 V 



(4.13) 
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4.2 Generalized correlation functions of SOS model 

For later convenience we also introduce another meromorphic function Qi n ) (w\Q kkl ''' k2n - :Lk ~ 21 , where 
^ I ^ n. Since k 2n — k — 21, the number of elements of the following set A' is equal to n: 

A' = AU {2n + 1, • • • ,2n + l}. 

Let k 2n +i = k 2n + 2i = k-2(l-i) for 1 < i < I. Then the meromorphic function ( w \Qkki-k 2n -ik-2l 
is defined as follows: 



' a-l 



2)1 



aeA 



Y[ { V a - U a + \ - k a } ( Y[ K ~ Uj - ^] ]J ^ ~ Va ~ ^ 

2n 

3=1 



1 

2 



2n+l r 1 1 7 1 

x -Q K »n+o -/■■„'} 

a'=2n+l 

2n-l 



[fa' - ^0 - 2J 

n^> _i n [««-«6+i]- x . 



(4.14) 



i<(j 



When I = 0, the meromorphic function (4.14) is evidently reduced to the expression of |UC| ). 
We further introduce the generalized correlation function of the cyclic SOS model as follows: 

dw a 

:,Jc a 2irv^lw a 



n 



kki ■■■k2n — ik—2l 



(4.15) 



where the integral kernel ty™ (w\C,) and the integral contour C a are respectively the same ones in the 
previous section. 

The formula Q4.1EQ can be used only for km ^ k. However, by noticing the fact that the Boltzmann 
weights of the cyclic SOS model ( 2.1C| ) is invariant under the shift (a,b,c,d) 1— > (— a. — b, — c, ~ d), we 
should obtain the expression for k2 n — k + 21 > k from (C)~ k ~ kl k 2n -i-k-2i ag f vj ows: 

f dWg 

In 2-Ky/^lWa 



■p( n ) ^Qkk 1 --k 2n -ik+2l 



n 

a&A'_ 



'i^^J l \w\QQ {n) (w\0 kkl - k2n - 1 



k+2l 



(4.16) 



where 

Q( n ) (yj\^\kkl — k2n-lk+2l 



1] R - U a + I + k a } [UK- Uj - ±] j[ [ 



It, - V a 



a£A- 



u' =1 



2n + Z r 1 1 1 7 1 

TT W - "0 + 2 + fc °'l 

/ 0.1 K' - u - 5] 

— 2n+l L 2 

2n-l 

nfer 1 n k-^+i]- 1 - 

j = l a,beA'_ 



Here 

A_ := {a\k a = fc _ a - 1, 1 < a < 2n}, A'_ = A_ U {2n + 1, • • ■ , 2n + /}, 
and A^n+i = ^2n — 2i = fc + 2(1 — i) for 1 ^ i ^ Z. 



(4.17) 
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4.3 First integral formula for the XYZ correlation functions 



Let 



op (Ci , • • • , c 2 «) = (v=i)- n e - w °) 

fe ,fel,---fc2n 



Here Cj = x Uj , and 



■(«an-tio)^(Ci 



X2n) kakl - k2n ■ 

(4.18) 



where F^\c) kkl '" k2n is denned by ( p^j4.14|) for k 2n < fc, otherwise by ( p^j4.17| ). The sum with 
respect to kj (1 ^ j ^ 2n) should be taken over fcj = ± 1. The sum with respect to kg should be 
as follows. When reQ and 2r = N/N' (N, N' are coprime), we have 



N+k 
N+k± 



1 (u) = (-iy +N 't k k±1 (u) 



is evidently invariant under the shift (k, ki, ■ ■ ■ , k2 n ) 



The cyclic SOS correlation function Fb- (C) 
(k + N, ki + N, ■ ■ ■ , + N). Thus the summand in ( 1.18] ) is also invariant under the same shift, so 
that the sum with respect to ko can be taken over Z/iVZ. Since there is no such invariance, the sum 
with respect to ko should be taken over Z when r is irrational. From these observations we have 



JV-l 

E if r e Q, 2r = N/N' (N, N' are coprime); 

= ( k =Q 
k & 



E 



(4.19) 



Simple observation shows that the i?-matrix symmetry (4.4) follows from the VF-symmetry (p. 2]) for 



( 4.15 4.16 ). The cyclicity for the generalized correlation functions in the cyclic SOS model does not hold 
but we have the following relations: 



Proposition 2 For fixed k = k, k\, ■ ■ ■ , k 2n _ 2 , k^n-i — with kj = kj-\ ±1 (1 ^ j ^ 2n — 1), the 
following cyclicity relations hold: 



s=±l 



skki ■■■k<2_ 



(4.20) 



s=±l 



[Proof] Let k' = k 
e = ± let 



21 + 1 with 1 ^ I ^ n. (The case k! — k + 21 — 1 can be similarly proved.) For 



h e (u) :-- 



'2- e f3+el27 



' u . 

v 2r ' 



(e<0) 
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Then from Q^9^ the relation fl4.2CQ is reduced to 



(-l) l [«2n -U + 1]J2 h e(k' ~ S(u 2n -U + 2))Fi")(C',X 2 C 2 „) 



kki ■■■k2n-2k k'+s 



s=±l 



[U 2n -U + 2]J2 sh e(k ~ S(u 2n - U + l))F^ (C,2n, C) 



k—skki ■■■kin- ik 



s=±l 



Here we used the abbreviation C' = (Ci> " ' ? Qin-x) an d the functional relation 



(4.21) 



/(u + 2) _ [tt+1] 



/(«) 



■2]" 



Simple calculations show that 



dw a 

,j Ca 1-K\T-Yw a 



s=±l 



n 



^ *(. ,l) (w ai ••• ,W 2rl+ l\(l,- ■ ■(2n)J r e,l(v\u ,Ui,- ■ ■ , U 2n ) , 



(4.22) 



where 



!F e ,t(v\uQ,ui, ■ ■■U2n) = h e (k -2l + 2 + u + u 2n - 2v 2n+ i) 

21 



[u 2n ~ Uq] TT { V 2n+i -UQ + \-k + 2(l~i)} 

TT 11 T ~ 11 [«2n+i - u 2r. 



[v 2n +l - "0 



2 J i=2 



[V2n+i —Uo — ^] 



' a-1 



2/i 



i] n 



«■> - v a 



j=a+l 



X {« a - U a + | - k a } Y[[v a -1 

■Ei \J = 1 

a<2n \ 

2n+l 2n-l 2n-2 

x (-i)'- 1 n n «i - u n^> _i n 

a'=2n+l j=l j = l a,6EA' 

a<i> 



(4.23) 



It follows from the definition that 2n £ A and A' = A U {2n +!,-•• , 2n + Z — 1} when s = 1, and that 



2n A and A' = AU {2n + 1, • • • , 2n + Z} when s = -1 in the LHS of (4.22). Thus we made the shift of 



variables (w 2n , ■ ■ ■ ,W2 n +l—i) (w 2n +i, • • • ,W2n+i) for s = 1 in (4.22). Furthermore, in order to derive 



the expression of (4.25) we used the identity 



=p- {v — mq — \ — k}h e (k + u — uq) — {v — u — ~ — k}h £ (k — u + uq) 



[V—Uo ,J 

[u — uo]{k}h E (k + 1 + uo + u — 2v) 
[v-uq- \] 



By repeating the similar calculations we find that eq. ( 4.21 ) is reduced to 



n 



- dl J ^ T — ^a L \ w \Ci, ■ ■ ■ C2n)Ai{v\u)B £ik .i{v 2n+ i,- ■ ■ , v 2n +l\uo, u 2n ) = 0, (4.24) 
,Jc a 2n^-lw a y 1 
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where 



a-l 



2n-l 



= Y\_ {v a ~ U a + ± - k a } \ [v a - u 2n - |] JJ K - u j - |] H 



a<2n 

2n+i 



3=1 



j=a+l 



n k 



wo 



i=l 
2n-2 



1 1 n ^ ~ ^ 2 " +i ~ 5] 



E.fc.i 



n^ -1 n K-vb+i}- 1 n [ Va - Va , - n K-^+ir\ 

j = l a,b£A> a,a'£A> a>,b>£A> 

a<b<2n a<2ji<a' 2n<a'<b' 

I 

h e (k -2l + u a + u 2n - 2v 2n +i) Y[{v 2n +i - Uo + § - k + 2(1 - i)}[v 2n+i - u 2n - ±] 

i=2 

l-l 

(-l) l h E (k -2 + u + u 2n - 2v 2n+ i) Y\_{v 2n +i - u Q + § - k + 2(1 - i)}[u 2n - v 2n+i - \}. 



i=l 



It is evident that (4.24) for I = 1 holds because B s ,fc,i = 0. In order to prove (4.24) for I > 1, we 
wish to show that 

i-i 

^e,fc,i( u l u 0,W2n) = ^ [l>2n+i ~ «2n+i+l + 1]B^ ; (w| U , U 2 n) , 



(4.25) 



where -B^, 1 is some function that is symmetric with respect to v 2n +i and v 2n +i+i- Suppose that (4.25) 
is true. Then AiB e ^.i is the sum of I — 1 terms, each of which is symmetric with respect to v 2n+ i and 
v 2n +i+i- Correspondingly, the integral (4.24) vanishes from the antisymmetry of with respect to 
u>'s. The claim of this Proposition is therefore reduced to eq. ( 4.25| ). 

Let us show ( [4.25 ) by the induction for I > 1. The validity of eq. (4.25) for I = 2 follows from the 
following identity: 



B £ ,k,2(v 2n +1, V 2n +2\u , U 2n ) 

= h e (k -4 + u + u 2n - 2v 2n+ i){v 2n+2 - u + § - k}[v 2n+2 - u 2n - |] 
+ h e (k - 2 + it + u 2n - 2v 2n+2 ){v 2n+ i - u + | - k}[u 2n - v 2n+1 - \] 
= -h s (k - 3 + u - u 2n ){k -3 + u + u 2n - v 2n+1 - v 2n+2 }[v 2n+1 - v 2n+2 + 1]. 



(4.26) 



For I > 2, the assumption of the induction implies 

1-2 

B ei fc-2,J-l(f2n+l! ' •• j V 2n +l-l\u Q , U 2n ) = ^2\v 2n+l - V 2n+l+ i + l]B^ k _ 2 j_ 1 (v 2n+1 , ■ ' ' , V 2n +l-l \uq, U 2n ). 



The validity of eq. (4.25) for general I thus follows from the relation 
B^k,i(v2n+i,--- ,v 2n+ i\u ,u 2n ) 

= B E} k-2,l-l(v2n+U ■ " j V 2n+ l-l\u 0l U 2n ){v 2n+ l - Uq + I - /c}[w 2 n+; " "2n ~ \] 

1-2 

+ (-l) l B e ,kA V 2n+l-l,V2n+l\ua,U2n) \\{V2n+i ~ Uq + § - k + 2(1 - i)}[u 2n - V 2n+i - ±], 



where we again use (4.26). □ 

As a Corollary of Proposition ||, we have 



Corollary 3 The cyclicity (£.1 ) holds for the XYZ correlation functions ({.lb). 
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4.4 Normalization in the XYZ case 



Furthermore, we find that the following normalization condition holds for the generalized correlation 
functions in the cyclic SOS model for fen = fe (mod 2): 



Fi" ) (C",C 2 „-i,^- 1 C 2 „-i) 



k"'k k ±\k% n 



= S k'M nT tj F -(C") k - k ' (* = ±), ( 4 - 27 ) 



where 



d, = 



e 2, 



(s > 0); 

~(l+r±2k) (s < Q) 



Proposition 4 The normalization condition j\4-(\ ) holds for the XYZ correlation functions ( 1.1$) . 
[Proof] First we notice the following identity 



' k-\~U . 7T^/— 1 \ £ f k — U . — 1 \ ft Q ( k — U . 7T V - 1 N 

7TV — 1 / U . 7rV— 1 \/j / k . TT\/~1 ^ 



C 2 [u){k}, 



where C is the same constant as the one defined in (£]§. Using ( jOf ) and (§J), we have 



J] tt.^-i-tloj^tit+'ftian- 



s=±l 



«2n=»2n-l-l 



-l{fc}5 £2ll _ 1+e2lli o- 



Furthermore, we also have 



s=±l 



— ,k 



4 +s (u 2n ~l ~ U f 2n -H k k +S {u2n ~ Uq) 



£2«e 2,. {k}S S2n _ 1+S2ni0 . 



(4.28) 



(4.29) 



(4.30) 



Thus, the normalization condition (4.6) in the XYZ Heisenberg antiferromagnet follows from (4.29) 



( pOD and ( p~27| ). □ 

From Corollary |^ and Proposition |4|, we have the following Theorem. 



Theorem 5 The XYZ correlation function Ci4"^(C) (4--^ ) solves the bootstrap equations 



4.5 Second integral formula for the XYZ correlation functions 



In this subsection we assume that r € Z>i. The integral formula (4.18) remains the sum with respect 



to the local height variables fe, k%, • ■ ■ , km in the cyclic SOS model. Thus, the expression (4. IS) consists 
of 2r x 2 2n terms because of (4.19) and the condition kj — kj-i ± 1 for 1 ^ j ^ In. 



Now we wish to present another integral solution to the bootstrap equations (4.4-4.6). Eq. ( 4.22 ) 
shows that the sum with respect to fen with fixed fe, fe, • • • , k^n—l can he performed to be in a simple 
form. In order to repeat this procedure with respect to fc2n-i) • • • fci, we find that the number of integral 
variables should be greater than In for the 2n-point correlation function. 
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Let kin = k + 2n — 21 (0^1^ 2n). Then the number of the set A is equal to 2n — I, where 

A := {a\k a = k a -i + 1, 1 < a < 2n}. 

We further introduce the set of indices 

A" = A U {2n + 1, • • • , (r + l)n + I}. 

Note that the number of the set A" is equal to (r + l)n. Let km+i — kin + 2i = k + 2(n — I + i) for 
1 ^ i ^ (r— l)n + l. Then the meromorphic function Q( n \w\£) kkl "' k2n with kin = k + 2n — 2l\s defined 
as follows: 



)(n) HC) /*... fc2 „ = JJ {Va _ Ua+ l_ ka} JJ[« - Uj --l] J] fo<-Vr-± 



o-l 2n 

*'j - "o - 2J 



(r+l)n+i r . x r 1 / 2n 



X 



2n-l 

n 1 n K-^+i]- 1 . 

j = l a,b£A" 

a<b 



Set 



F(n) {c) k kl - k2n = ~ n Y[ ^g{ Zj /z k )xF^\o kk -- k ^, (4.32) 



where c„ is some constant, and 



•„ J Ca 2TTy/^lw a 



Here, the kernel has the form 



2n 

a"G/li=l J lsjj<fc^2n 

The function ^"''(HC) is a function of (r + l)n w's and 2n £'s. The properties of ^i™** (HO is the same 



as those of 1?^ i but in the RHS of ( 3.14 ) the product with respect to a should be taken over A" . 



The second integral formula of the XYZ correlation functions is as follows: 

Gi n \() = (V=T)-" Yl *fe(«i-«o)®---®*t;- l (u an -uo)FW(C)* * 1 "-* a -. (4.35) 



k Q, k l,--- k 2n 



This formula essentially solves the first two of the bootstrap equations. The i?-matrix symmetry (4.4) 
evidently holds. The modified cyclicty 

P12 • ■ •P2„-i2„G^ ) (C'^ 2 C2„) - <r(-l) (r - 1)n G£°(C»n>C), ( 4 - 36 ) 
can be proved in a similar manner as in Proposition @ and Corollary a. 
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The advantage of the second formula consists in the fact that the sum with respect to ki , • • • , k 2n 



can be carried out. By repeating the similar calculation as (4.22) was derived, we find 

dw a 



Gi n Ho= n / 

aeA" JC 



where 



2r-l 



QW(ttlO Ei '" E2n = ( A /^i)-"+£ 3 =o £ j( fc ~ 2 »^ i +3)(-i)( r + 1 )" 

c i / a—1 2n 

K ~% iik-^-2-] n b 



X 



fc=0 
2n 

n 



(r— l)n r i , _ , 

n{«2«+j - u + o - fc - 2(n + z)} 




2 



J 2 J 11 1"J 

j=a+l 

2n 



V a ~ h] 



(4.37) 



( = 1 
2d 



[V2r. 



Uo 



n 



V 2n +i - Uj - 2J 



Y\_f(uj -u )h ej (k+j + u + Uj -2vj) [u a - u 6 + 1] 



»,b£ A'- 



Finally, we have to carry out the sum with respect to A; of the following form: 

2r-l 2n 

F^(u \{u J } ls : J<2n \{v a } 1 ^ {r+1)n y = ^2{V^l) k ^=° e3 Y[h ej {k+j + u + u j -2v j ) 



fc=0 
(r-l)r 



J| {«2n+ l - Uo + I - k - 2(n + ?:)}. 



(4.38) 



This function has precise quasi-periodicities as a function of Uq, Uj (1 ^ j ^ 2n) and w a (1 ^ a ^ 
(7' + l)n). In particular, you can easily find the Mo-dependent part of F^. Furthermore, if the limit 
r — > 1 is taken, the explicit expression for n = 1 can be obtained as follows: 

F^(u \u 1 ,u 2 \v 1 ,V2) ElE2 

C _e 2 ^^-v 1 +v 2 i^)e 3 (u + ^±^-v 1 -v 2 -,^) (£i£ 2 <0); (4.39) 

C+B x (21=31 - Wl + W 2 ; SSO ) 4 (tio + Hi±Hl _ Vl _ „ 2 . SO ( £l£2 > 0). 



Here C± are some constants times exponential functions of uq, u\, u 2 , Vi, v 2 . 



5 Concluding remarks 

In this paper we have constructed two integral solutions to the bootstrap equations for the XYZ Heisen- 
berg antiferromagnet. These solutions are expected to give the correlation functions of the XYZ model 
in the antiferromagnetic regime. The first solution is essentially the same as the formula given by 
Lashkevich and Pugai J^] . Lashkevich-Pugai's formula can be obtained from the correlation function in 
the RSOS type model, by using the vertex-face correspondence. In order to avoid the pole resulting from 
[k] in the denominator, k £ Z + 5 should be assumed with some real <5, and the limit (5^0 should be 
taken after all calculation P] . On the other hand, we constructed our first formula from the correlation 
function in the cyclic SOS model so that no such regularization was needed. 
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By the construction, both Lashkevich-Pugai's formula and our first formula for 2n-point XYZ cor- 
relation function remain the sum with respect to the local height variables ko , k\ , • • • , k^n other than 
n- fold integral. When r G Z>i we construct another (r + l)?i-fold integral formula, in which the sum 
with respect to k±,- • • , kin can be carried out. The structure of our second formula is quite similar to 



Concerning the second formula, we have not yet proved the normalization condition. The number of 
integral variables in the second formula are much greater than that in the first one, so that the recursion 
relation for the second one is not so simple. Nevertheless, we believe that the normalization condition 
holds in the second integral formula. Actually, the number of integral variables in the second formula 
is not minimal but redundant. It is sufficient for the number of integral variables N in the second 
formula to satisfy n ^ N — n = (mod r). Thus, we can make an A r (n)-fold integral formula, where 
N(n) = n + mr for (m — l)r < n ^ mr with some m £ Z >0 . If we employ this 'minimal' integral 
solution, the normalization condition holds unless n = 1 (mod r) from the same reason as in the first 
formula. However, the validity of the recursion relation in the 'minimal' integral solution becomes again 
unclear for n = mr + 1 because N(n) — N(n — 1) = r + 1. 

Another future problem is to consider the bootstrap equations for Belavin's Z/nZ-symmetric model 
and to construct integral formulae of the correlation functions. We wish to address this problem in a 
separate paper. 
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A Properties of the i?-matrix 

The main properties of the i?-matrix of the eight-vertex model are the Yang-Baxter equation 



where the subscript of the R- matrix denotes the spaces on which R nontrivially acts; the initial condition 



Shiraishi's formula p3j, though the latter treated only the case r 



3 
2 ' 



J Rl2(Cl/C2)^13(Cl/C3)i?23(C2/C 3 )=i?2 3 (C2/C3)i?13(Cl/C3)i?12(Cl/C2); 



(A.l) 



R(1) = P; 



(A.2) 



the unitarity relation 



i?i2(Ci/C2)i?2i(C 2 /Ci) = i; 



(A.3) 



the Z2-parity 



fll2(-C) 



(A.4) 
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and the crossing symmetries 

J?2 1 i(C 2 /Ci) = <7^i2(^- 1 Ci/C 2 )<7f = -(TyJZxaC-z-^x/CaK. (A.5) 



In (A. 4, A. 5) the shift £ l_ > — C implies the one such that u i— ► u — ^ — . The properties (A.2-A.5) hold 



if the normalization factor of the R- matrix satisfies the following relations: 

K(o«(r x ) = i. ^r 1 ) = «(0- (a.6) 

Under this normalization the partition function per lattice site is equal to unity in the thermodynamic 
limit @ ||. 
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